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In scale invariant hydrostatic barotropes, the radial evolutionary equation linearly relates the local 
gravitational and internal energies. From this first-order equation, directly follow all the properties 
of polytropes and the important mass-radius relation. Quadrature then leads to the regular Lane- 
Emden functions and their Picard and Pade approximations, which are useful wherever stars are 
approximately or exactly polytropic. We illustrate this particularly for the n=3 regular polytrope 
and obtain analytic approximations to the solution of the Lane-Emden equation, valid over the bulk 
of relativistic degenerate stars (massive white dwarfs) and chemically homogeneous stars in radiative 
. equilibrium (ZAMS stars). 
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I. DIFFERENT INDEX POLYTROPES SHOW DIFFERENT CORE CONCENTRATIONS & 

MASS-RADIUS RELATIONS 



In the preceding paper [1], we derived the characteristic differential equations 

d\ogv n dlogu 

7— — = o =dlogr 1 

u — 1 + v n 3 — u — nv n 

for the homology invariants u := d log m/d log r, v n := — dlogP/ p/dlogr of n-th order polytropes. Each regular 
polytrope is a solution v n (u) of the first-order Abel equation, subject to the regularity condition u(0) = 3,v n (0) = 

dv n Vn(u- ,<y, 

du u(3 — u — nv n ) 
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characterized by the boundary value oU n := cj n (£i n ) of its homology invariant 

"»K):=K)* /2 • (3) 

As tabulated in the third, fourth and last columns of Table I 0-01? o^n determines the boundary £i n and the density 
ratio pcn/pni where p n = 3M/4:7rR 3 is the mean density, the M-R relation M ~ QU n 

R (n-3)/(n-i) m These properties 

of polytropes result directly from the first-order equation (2), which encapsulates the effect of scale invariance, and 
whose solution for each polytropic index n is plotted in Figures 1 and 2. 



A. Regular Polytropes of Different Index Have the Same Inner Core But Different Envelopes 

After solving the second-order Lane-Emden equation, each Lane-Emden function of index n is also characterized 
by the first zero of n (£)> the dimensionless radius £i n . The sixth and seventh columns in Table I list dimensionless 
values for the inner core radius r- lcn /R = ^ cn and included mass mi cn /M, shown by red dots in Figures 3, 4. The 
inner core radius £i cn , defined implicitly by u(^ cn ) := 2, is where the acceleration Gm/r 2 reaches a maximum and the 
gravitational energy density overtakes the internal energy density. In the inner core, the internal energy dominates; 
in the envelope outside, the gravitational energy dominates. 

For homology variables, we prefer a new independent variable z := 3 — u = —dlogp n /dlogr and a new dependent 
variable w n (z) := nv n := — dlog p/dlogr. In term of these invariants, the characteristic differential equations (1) are 

dz d\ogw n 

— = diogr . (4) 



(3 - z)(w n - z) 2 - z + 

The first equality is the first-order equation (2) for the invariant w n {z), which we solve for the central boundary 
condition w n — > 5z/3 when z —> 0. w n (z) and the differences (3/5)[w n (z) — w^(z)] = (3/5)w n (z) — z = dlog p/ps ~ G n 
are plotted in Figures 2 and 1, respectively, for polytropic indices n — 1,2, 3, 4, 5. 

For incompressible matter (n = 0), there is no core concentration: the mass is uniformly distributed, and the entire 
star is core. But as the equation of state softens as n increases toward 5, the gradient decreases, the core concentrates, 
the inner core radius shrinks, and the envelope outside the core grows: r- lcn / R — >• 0, mi cn /M — 0.19. For the softest 
equations of state n < 5, the stellar radius £i n « 3(n + 1)/(5 — n), the inner core radius shrinks £i cn w yl0/3n, their 
ratio r icn /R = 6cn/6n « 0.045(5 - n), m icn /M w 0.20, and uj n w \f3/^ ln . 

For n < 5, the stellar boundary lies at finite radius. The Noether charge is nearly conserved at G n ~ in the inner 
core, but grows rapidly as the boundary is approached (Figure 1). The polytropic form is locally scale invariant, but 
the radius R determines the stellar scale. 

For n = 5, the core becomes infinitely concentrated, shrinking to zero, and the star is all envelope. The n = 5 
regular solution 

e 5 (o = (i+em- l/2 , p = Pc (i+em- 5/2 , m = M-e/(z+e? 12 , « 6 = (3-«)/3 . (5) 

= 0, so that the n = 5 poly trope has infinite stellar radius R for any mass M and is globally scale invariant. The 
Noether charge G5 ~ (1*5 + u/3 — 1) =0 everywhere. 



TABLE I: Scaling Exponents, Core Parameters, Surface Parameters, and Mass-Radius Relations for Polytropic Gas Spheres of 
Increasing Core Concentration 
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FIG. 1: Dilution of polytrope density profiles in the envelope as the boundary is approached (u — Y 0). All solutions approach 
the same density structure w n (z) — »• 1^5 = (5/3) (3 — it) at the center (it — >> 3), but differ outside the core. As the boundary is 
approached, the mass m — » M, and w n n[ocj^ _1 /^] 1//n , r/iUn — > 0. 



B. Lane-Emden Functions and Their Picard Approximations 



For < n < 5, 



M = f 
Jo 



dz w n f Z -± Wn,n) « (5/J)[l - (1 - z/i) J ] := W » Pic (z) , J:=(9n-10)/(7-n) , (6) 



is well- approximated by the Picard approximation, obtained by inserting the core values w n (z) « (5/3)z inside the 
integrals. Indeed, this Picard approximation is everywhere exact for n = 0, 5. For < n < 5, it breaks down 
only in the outer envelope, where w n diverges as w n — >> n^cj^ -1 /^] 1 / 72 . Integrating over z, the density profile and 
Lane-Emden functions are [7[ 



Pn(z)/Pcn = K = ex P 

# n = exp { - / — — 

^ JO n l w n 



f_ r dz Wn (z) i 5/2 



(z)-z](3-z). 



(8) 
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FIG. 2: Different Lane-Emden density gradients w n (u) := — d\ogp/d\ogr as the outer boundary u = is approached. Inside 
the core u < 2, all w n ~ if 5 = (5/3) (3 — u). For the softest equation of state n — 5, the stellar radius is infinite, the density 
gradient w*>(u) — (5/3) (3 — u) everywhere, and p — p^ 5 p 3 ^ 5 . For stiffer equations of state n < 5, the density gradients w n (u) 
steepen as the finite radius is approached u —> 0. 



where again the Picard approximations are obtained by inserting the core relations w n (z) = nv n (z) ~ (5/3)z inside 
the integrals. 



II. IN P — Kp 4/S POLYTROPES, RADIUS IS INDEPENDENT OF MASS 

The n = 3 polytrope, which is realized in relativistic degenerate stars (massive cold white dwarfs) and in the 
Eddington standard model (luminous ZAMS stars), is distinguished by a unique M — R relation: the mass M = 
47r 0^3 (K/ttG) 3 / 2 depends only on the equation of state constant K := P/p 4 / 3 , where the constant K depends on 
the entropy but is independent of the radius R. In these stars, the gravitational and internal energies cancel, making 
the total energy W = ft + U = 0. Such stars are in marginal dynamical equilibrium at any radius. 
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FIG. 3: Normalized density profiles as a function of fractional included mass m/M, for poly tropes of finite mass M and 
compressibility increasing with n. The red dots mark the cores. For incompressible matter (n = 0), the polytrope is all core. 
As the matter softens (n increases), an envelope grows to ultimately encompass just over 80% of the mass. For any n > 1, the 
density at the inner core radius stays in the narrow range 0.37 < p(n C n) / pc < 0.42. 



A. Relativist ic Degenerate Stars: K Fixed By Fundamental Constants 



Massive white dwarfs are supported by the degeneracy pressure of relativistic electrons, with number density 
n e = p/fi e mH, where ra# is the atomic mass unit and the number of electrons per atom p e = Z/A = 2, because 
massive white dwarfs are composed of pure He or C 12 /0 16 mixtures. Consequently, Kwd = [f Y^ m HPe~^^ 
depends only on fundamental constants. This universal value of Kwd leads to the limiting Chandrasekhar mass 
M Ch = ^MJ£ = 5.83M / M 2 = i.46M . 



B. Zero-Age Main Sequence Stars: K(M) Depends on Specific Radiation Entropy 

In a non-degenerate ideal gas supported by both gas pressure P gas = 1ZpT / p and radiation pressure P ra d = aT 4 /3, 
the total pressure P = P gas + P rad = P gas [l + (1 - ft)/ ft] = P gas [l + pS rad /47l], where /3 := P gas /P and 

4r/T 3 

S™d = — = (ft//*) • 4(1-/3)//? , 5 gas = (K/vl) ■ log (T 3 / 2 /p) (9) 

are the specific radiation entropy and gas entropy for an ideal monatomic gas. Bound in an n = 3 polytrope, the 
gas entropy is a constant depending on K, the total specific entropy is constant, and the total energy W = O + U 
vanishes. Because such stars are in marginal dynamical equilibrium at any radius, their mass M is independent of 
radius R. 

In zero-age Main Sequence (ZAMS) stars of mass M, the radiation entropy 5 ra d is approximately constant, because 
at each radius, radiation transport leaves the luminosity generated by interior nuclear energy generation proportional 
to the local transparency (inverse opacity) n~ x . Eddington's standard model assumes constant S ra d(M). This makes 
/3(M) and K(M) = {[3(1 — (3) /a](7Z/ pfS) 4 "} 1 / 3 constants depending only on M, according to Eddington's quartic 
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FIG. 4: Normalized density profiles as function of fractional radius r/R. The density is constant for incompressible matter 
(n = 0) which is all core, but is concentrated at the origin for an unbounded poly trope (n = 5, R = oo) which is all envelope. 
For any n > 1, the density at the inner core radius stays in the narrow range 0.37 < p(r\ cn )/ p c < 0.42. 



TABLE II: Taylor Series and Picard Approximations n pic to Lane-Emden Functions n 



n Lane-Emden Function and Taylor Series 



N := 5/(3n - 5) Picard Approximation 6> n pic := (1 + £ 2 /6iV)" 



1 - £ 2 /6 

sinf/f - 1 - £ 2 /6 + £ 4 /120 - f 6 /5040 + • • • 
1 - £ 2 /6 + n£ 4 /120 - n(8n - 5)/15120£ 6 + • 

a+e 2 /3)- 1/2 



-i 

-5/2 

5/(3n - 5) 
1/2 



1 - f/6 

(l-e 2 /l5) 5/2 : 

(1 + ^/6^)-^ 
(l + £ 2 /3)- 1/2 



1 - £ 2 /6 + £ 4 /120 - ^ 6 /!0800 + • • • 
= 1 - ^ 2 /6 + < 4 /120 - n(6n - 5)£ 6 /10800 + • 



equation 



, 2 x2 



1^ - ) A I 



3VT0 CJ3 / ^ \3/2 

^ 3 V^7^J 



18.3M 







(10) 



The luminosity L = LEdd[^ — P(M)] depends on the Eddington luminosity LEdd '•= ^kcGM / k v through the 
photospheric opacity k p . From Eddington's quartic formula, the stellar luminosity 



L/Lq = ^^ {0 .oo3)^\M/M Q f. 



(ii) 



This is confirmed [5] in ZAMS stars, particularly over the lower end of the Main Sequence O.IMq < M < IOOMq. 

III. ANALYTIC APPROXIMATIONS TO LANE-EMDEN FUNCTIONS 

We now find analytic approximations to the radial structure of any poly trope, particularly the important n = 3 
polytrope. From the solutions w n (z) to the first-order equation, we now use 



dm/m := u • dr/r = dz/[w n (z) — z] 



(12) 
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to obtain 



„,W/M = ,!^.e X p{/ o d5 {_L_ <£)»/* (13) 

rUVR - (|)" 2 ■ e XP {/' ^{ (3 _ a)[ J n(z) _ zl 4}} » ^ (14) 



for the mass and radial distributions. The integration constants R, M, p c express the scale dependence of the 
polytrope. Using equation ([T2]) to eliminate the Picard approximations 

^nPic(0 = (l + ^ 2 /6iV)- iV , 7V:=5/(3n-5) (15) 

to the Lane-Emden functions are obtained and tabulated in the last column of Table I. For n = and 5 polytropes, 
this closed form is exact. For intermediate poly tropic indices < n < 5, the Picard approximation breaks down 
near the outer boundary, but remains a good approximation over most of the polytrope's bulk. The worst Picard 
approximation is for n = 3. 

For n = 3, the Taylor series expansion 

#3(0 = 1 - 0/6 + 0740 - (19/5040)0 + (619/1088640)0 - (2743/39916800)0° + • • • (16) 
= 1 - 0.16666670 + 0.0250 - 0.00376980 + 0.00056860 - 0.000068720° + • • • , (17) 

when truncated at tenth-order, has a radius of convergence at £ ~ 2.4. Outside this radius, the Picard approximation 
(15) 

03Pic(O = 3 +0/3528 + ... (18) 

overestimates the Taylor series expansion, but remains a good approximation over the bulk of the star, with error 
< 10% out to £ w 3.9, more than twice the core radius and more than half-way out to the stellar boundary at 
£13 = 6.897. The Picard approximations in white dwarf and ZAMS stars suffices, except for their very outer envelopes, 
which contain little mass and are never poly tropic. 

Because it satisfies the central boundary condition but not the outer boundary condition, the Picard approximation 
underestimates 0'(£) and overestimates 0(£) outside £ ~ 3.9. A much better and simpler approximation is the Pad 
rational approximation Q 

3 Pad := i JJ^i^alcn.^o = l--0.166667e 2 + 0.025e 4 -0.00376984e 6 + 0.0005686e 8 -0.0000857618e 10 +---, 
1 + 119^/756 + £ 4 /10(J8 

(19) 

which, already in tenth-order, starts converging faster than the Taylor series. In fact, this Pade approximation agrees 
almost exactly with the exact solution out to £i = 6.921, very close to £13 = 6.897, the first zero of the exact Lane- 
Emden function 0${£). These simple analytic approximations to 03(£), shown in Figure 5, simplify structural modeling 
of massive white dwarfs and ZAMS stars. 

IV. CONCLUSIONS 

For spherical hydrostatic systems obeying barotropic equations of state, the scaling evolutionary equation leads 
to a linear relation between the local gravitational and internal energies. From this equation, directly follow all the 
properties of polytropes and the important mass-radius relation. Quadrature then leads to the regular Lane-Emden 
functions and their Picard and Pad approximations, which are useful wherever stars are approximately or exactly 
polytropic. 

We illustrated this for the n = 3 regular polytrope and obtained analytic approximations to the solution of the 
Lane-Emden equation, valid over the bulk of massive white dwarfs and ZAMS stars in radiative equilibrium. 
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FIG. 5: The exact Eddington Standard Model Lane-Emden function 03(0 and its Picard and Pade approximations. Even in 
this worst case, the Picard approximation works well out to twice the core radius at £i C 3 = £i3(nc3/-R) = 1-65, but breaks 
down near the boundary. The Pade approximation is indistinguishable from the exact solution, vanishing at £i = 6.921, very 
close to the exact zero £13 = 6.897 (For polytropes, the Picard approximations would be even better and become exact 

everywhere as n — >> or 5.) . 
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